arXiv:1503.01710vl [physics.space-ph] 5 Mar 2015 


JOURNAL OF GEOPHYSICAL RESEARCH, VOL. ???, XXXX, DOL10.1002/, 


Resonance broadening due to particle scattering and 
mode-coupling in the quasi-linear relaxation of 
electron beams 

1 1 .12 
Nicolas H. Bian , Eduard P. Kontar , and Heather Ratcliffe ’ 


Corresponding author: N.H. Bian School of Physics and Astronomy, University of Glasgow, 
G12 8QQ, Scotland, UK (n.bian@physics.gla.ac.uk) 

X SUPA, School of Physics and Astronomy, 

University of Glasgow, G12 8QQ, Scotland, 

UK 

2 Centre for Fusion, Space and 
Astrophysics, Department of Physics, 

University of Warwick, CV4 7AL, UK 


DRAFT 


March 6, 2015, 2:35am 


DRAFT 



X - 2 


BIAN ET AL.: RESONANCE BROADENING IN QUASI-LINEAR RELAXATION 


Abstract. Of particular interest for radio and hard X-ray diagnostics of 
accelerated electrons during solar flares is the understanding of the basic non¬ 
linear mechanisms regulating the relaxation of electron beams propagating 
in turbulent plasmas. In this work, it is shown that in addition to scatter¬ 
ing of beam electrons, scattering of the beam-generated Langmuir waves via 
for instance mode-coupling, can also result in broadening of the wave-particle 
resonance. We obtain a resonance-broadened version of weak-turbulence the¬ 
ory with mode-coupling to ion-sound modes. Resonance broadening is pre¬ 
sented here as a unified framework which can quantitatively account for the 
reduction and possible suppression of the beam instability clue to background 
scattering of the beam electrons themselves or clue to scattering of the beam- 
generated Langmuir waves in fluctuating plasmas. Resonance broadening be¬ 
ing essentially equivalent to smoothing of the electron phase-space distribu¬ 
tion, it is used to construct an intuitive physical picture for the stability of 
inverted populations of fast electrons that are commonly observed in-situ to 
propagate in the solar-wind. 
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1. Introduction 

Solar flares are transient bursts in the solar atmosphere following the release of energy 
accumulated in the coronal magnetic held. A major aspect of solar hares is the high 
acceleration efficiency of electrons inferred from their hard X-ray (HXR) emission [see e.g. 
Holman et al, 2011; Kontar et al, 2011a, for recent reviews] . Although hare accelerated 
electrons are prominently diagnosed via their HXR emission at the dense footpoints of 
coronal magnetic loops, these electrons can also escape into interplanetary space and be 
detected in-situ near-Earth or remotely via their radio emission in the solar wind [e.g. 
Lin, 2011, as a recent review]. As a population of escaping electrons travel further away 
from the Sun along the interplanetary magnetic held, faster electrons outpace the slower 
ones creating an electron beam with positive gradient in velocity space. When the number 
of beam electrons is sufficiently large, the beam becomes unstable to the generation of 
Langmuir waves which then convert to electromagnetic waves [ Ginzburg and Zhelezniakov , 
1958]. This plasma emission process can be seen via type III radio bursts. The rapid 
frequency drift of these bursts is associated with the beam travelling away from the Sun 
in a plasma of decreasing density, see e.g. Goldman [1983]; Melrose [1985]; Robinson and 
Cairns [1998]; Melrose [2009] for reviews. A fundamental problem in the theory of type 
III solar radio emission is therefore the understanding of the interaction of the exciter 
with the background plasma which is not quiescent nor homogeneous. For instance, 
turbulent magnetic fluctuations in the solar wind are known to efficiently scatter non- 
thermal particles resulting in broad pitch-angle distribution despite the effect of adiabatic 
focusing of the held, see e.g. Shalchi [2009] for a review. These turbulent magnetic 
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fluctuations also result in cross-field transport of flare accelerated electrons in the corona 
[Bian et al, 2011; Kontar et al, 2011b] and in the solar wind [ Shalchi, 2009]. Beam 
generated Langmuir waves can also be efficiently scattered by density inhomogeneities. 
The effect of turbulent density fluctuations has been widely discussed in the context of 
the beam-plasma relaxation. Vedenov et al. [1967] first described refraction by random 
large-scale density inhomogeneities as a diffusive transfer of Langmuir wave energy in 
wave-number space. This diffusion equation in wave-number space was then used to 
describe the spectral transfer of wave-energy out of the region of excitation by the beam 
[Nishikawa and Ryutov, 1976; Muschietti et al, 1985] providing a path for the weakening 
of the beam instability and also for acceleration of the electrons during beam-plasma 
relaxation [ Escande, 1979; Kontar et al, 2012; Ratcliffe et al, 2012], Beam-generated 
Langmuir waves are shifted toward higher/lower phase velocities when propagating into 
a region of higher/lower plasma density so that electrons with higher/lower velocities 
can interact with these waves [Breizman and Ryutov, 1969]. Mode-coupling of Langmuir 
waves, for instance to sound waves, can also play a stabilizing role on the beam plasma 
system, Papadopoulos et al. [1974]; Rowland and Papadopoulos [1977] have suggested that 
the modulational instability can remove beam-generated Langmuir waves from resonance 
with the beam by scattering them to different wavenumbers. The important role of 
density fluctuations is attested by the clumpy spatial distribution of the wave field which 
is measured in situ in the solar wind [ Gurnett and Anderson, 1977; Lin et al, 1981] as 
such localization is probably associated with density cavities [see e.g. Zaslavsky et al, 
2010; Hess et al, 2010]. 
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In this work we consider anew this important problem of beam relaxation in turbulent 
plasmas. We show that resonance broadening provides a unified framework to account for 
the effects of turbulent scattering of the particles and scattering of the waves in the quasi- 
linear relaxation of electron beams. Resonance broadening was first discussed by Dupree 
[1966] in the context of quasilinear theory in order to treat certain features associated 
with the high-intensity of the wave electric field, see also Rudakov and Tsytovich [1971]. 
Independently of its origin, resonance broadening filters out the small scale variations of 
the electron phase-space distribution and therefore, it provides an intuitive picture for 
the weakening and possible suppression of the beam instability due to particle and wave 
scattering. 

We start by reviewing in Section 2 major aspects of the quasilinear description of re¬ 
laxation of electron beams. In Sections 3 and 4, it is shown that the effect of scattering 
of the particles and waves can be included in the form of broadening of the wave-particle 
resonance involved in the quasilinear diffusion equations. In Section 5, we start from the 
Zakharov [1972] equations and derive a resonance broadened version of weak-turbulcnce 
theory, where broadening of the wave-particle resonance occurs as a result of non-local 
three-wave coupling between two Langmuir (L) modes and an ion-sound (S) mode with 
^l, kjj 3> ks- A summary of the results and conclusions are given in Section 6. 

2. Quasilinear diffusion equations for resonant wave-particle interactions 

Quasilinear theory is one of the few theories of plasma turbulence [Sagdeev and Galeev, 
1969; Davidson, 1972; Melrose, 1980; Tsytovich, 1995]. It is central to the modelling of 
stochastic acceleration of particles by waves and plasma turbulence, for instance during 
solar flares [see e.g. Miller et al, 1997; Petrosian, 2012; Bian et al, 2012, as reviews]. 
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Originally, quasilinear theory was proposed as a mean-field description of Langmuir wave 
generation by a population of beam-electrons, also accounting for the self-consistent evo¬ 
lution of the wave spectrum that these electrons may emit and reabsorb in the course 
of the bump-on-tail instability [Pines and Schrieffer, 1962; Vedenov and Velikhov , 1963; 
Drummond and Pines, 1964], The quasilinear diffusion equations describing the relax¬ 
ation of beam electrons are commonly derived from the Vlasov-Maxwell equations, which 
for a one dimensional electron plasma read 

m +v m + ^ {x ' t) iL] f{x ’ v ’ t) = 0 ’ (1) 

i>i: . 'I f = -47T en j dvvf{x, v, t). (2) 

Here, f(x,v,t ) is the electron distribution function (normalized such that / dvf(x,v,t ) = 
1), E(x,t ) is the electric held, and e, m e are the electron charge and mass respectively. 
In this high frequency limit, the ions are forming a static neutralizing background and, 
hence, the electric current j = e f dvvf(x,v,t ) is carried entirely by the electrons. 

Let us start by considering the problem of the statistical acceleration or deceleration of 
electrons in a broad spectrum of randomly varying longitudinal electric fields [ Sturrock , 
1966]. The Newton equations of motion for the electrons under the influence of electric 
fields are 


dv 

dt 


— E (x,t ) , 
m e 



( 3 ) 


The electric held is assumed to be statistically homogeneous and stationary, with the 
property (E) = 0 and (E 2 ) ^ 0. For the stochastic acceleration problem at hand, the 
fluctuating electric held is also characterized by a a correlation length A and a correlation 
time r. The evolution of the particle distribution function may then be given by a Fokker- 
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Planck equation, 


df(v,t) = d df(v,t) 

dt dv dr ’ U 

where D is the coefficient involved in the diffusion of electrons in velocity space. Under 
this diffusion, particles flow down the gradient in velocity space. This tends to flatten this 
gradient of distribution function and, hence, the action of randomly varying electric fields 
leads to statistical acceleration of the particles when df /dv < 0, or to their statistical 
deceleration when df /dv > 0. Accordingly, there is damping or growth of electric energy. 
The quasilinear approximation is a perturbation scheme based on the smallness of the 
so-called Kubo number, which is essentially the normalized amplitude of the electric field. 
Such perturbation scheme is used for obtaining an expression for the diffusion coefficient 
D entering Eq.(4) in terms of the spectrum of the electric field fluctuations. The self- 
consistent rate of damping or growth of the electric field is then determined by applying 
conservation of total energy, kinetic plus electric in the system, in the particular case of 
waves having specific dispersion relation ov(k) relating the frequency and wave-number 
involved in the definition of the Fourier components of the electric held fluctuations. This 
proceeds as follow. The Fourier components E k ^ of E(x,t ) are defined via 

(5) 

k U 

with k = nSk, u = m5uj, 5k = 27 t/L, and 5lu = 2n/T. ffere, n and m are integers while L 
and T are the spatial and temporal periods of the electric held. Under the Corrsin [1959] 
independence hypothesis, the diffusion coefficient in velocity space can be expressed as 
[Krommes , 2002] 


D 


mi 


XXI AT 

k u 



dt{e ikx ~ iut ) . 


( 6 ) 
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In the limit of densely packed Fourier harmonic, one is lead to replace the discrete sum 

by a continuous one Yhk,u J dk dui/Sk Sui , and since by definition the electric field 

spectrum is given by 

S E (k,u) = \E k>u \ 2 /6kSu, (7) 

one also obtains that for a continuous spectrum of fluctuations, 

D = —l [ dkduS E (k,uj) [°° dt (e ikx ~ iu}t ) . ( 8 ) 

mj J J Jo 

In this expression for the diffusion coefficient in velocity space, the function given by 

roo 

R(u-kv,Auj)= dt(e ikx - iut ), (9) 

Jo 

is the resonance function which characterizes the wave-particle interactions. A more 
precise meaning of the averaging procedure (...) and of the resonance width A ui involved 
in this expression for the resonance function will be given below. 

We notice that the equations of motion for the particles can also be written in the 
following non-dimensional form : dv/dt = eE, dx/dt = v. Time and length have been 
normalized to r and A and E has been normalized to yj(E 2 ). The Kubo number e = 
er 2 ^1 (E 2 ) /m e A represents the normalized amplitude of the electric field, so that when e <C 
1, the resonance function (9) can be calculated via the unperturbed particle trajectory. 
This unperturbed trajectory is given by x = vt and v = const which are the solutions of 
the equations of motion for e = 0, i.e. in absence of the field. The small Kubo number 
limit corresponds to the weak-held limit discussed by Sturrock [1966]. In absence of 
background scattering, the e = 0 zeroth-order particle trajectory is deterministic, hence 
we also drop the average involved in the resonance function (9). In this resonant case, the 
integral given by f^°dte^ kvt ~ ujt ' > is not convergent when t —> oo, so that an infinitesimal 
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damping factor v = 0+ is added in the exponential, 

dt e i(kvt-ut) e -vt = -* - M 7 T 6(u - kv) + iV ( —) , (10) 

(kv - uj) + iv y ’ \u-kvj y ’ 

where V is the principal value. Only the real part contributes to the resonant diffusion 
coefficient (8) which takes the form 

D = —— f j dkdu SE(k,uj) 5(lv — kv). (11) 

Tfl e J J 

A Fourier component [tu, k\ of the electric field fluctuations can exchange energy only with 
particles having velocity v such that the resonance condition v = uj/k is fulfilled. 

The only linear modes in Equations (l)-(2) are Langmuir waves. Therefore, we now 
specialize to electric field fluctuations which are produced by Langmuir waves whose 
dispersion relation is u(k) ~ u pe = ^/47rn e e 2 /m e , and write SE(k,ui ) = Ssik) 5(lu — u pe ). 
It is not here important to include the thermal corrections to the dispersion relation 
for k\oe 1. The spectral energy density of Langmuir waves is given by W(k) = 
2(S l £:(/c)/87r), where the factor of two accounts for the electric energy and the kinetic 
energy of the thermal electrons participating in the wave motion. This leads to the 
following diffusion equation 



df(v,t) 47r 2 e 2 d W(k — u pe /v,t) df(v,t) 


dt 


ml dv 


dv 


( 12 ) 


describing stochastic acceleration of electrons by a spectrum of Langmuir waves. Let us 
come back to the case of a single Fourier component in the electric field E k ^ e A kx ~ ut ) +cc, 
where cc denotes complex conjugate, so that the diffusion equation is 


<9/CM) 7re2(9 |A |2 r/ , , d f( v ,t) 

\E k ^\ 6(u-kv)- 


dt 


ml dv 


dv 


(13) 


For such a single Fourier component, the sharp resonance function <5(ca — kv) involved in 

the diffusion equation (13) implies diffusion of the electron distribution function f(v) over 
DRAFT March 6, 2015, 2:35am DRAFT 
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a single point in velocity space where v — u/k. This is a singular and ill-defined diffusion 

process. A main consequence of the resonance broadening effect to be discussed below, 
is that when the resonance function has a finite width Ac o, even a plane wave can now 
produce diffusion over a broad region in velocity space. This region in velocity space is 
centered around v = to/k and has an extent Av = Au/k where Acu is the broadening 
width. However, we know that a plane wave undergoes damping while accelerating the 
particles, and therefore, the singular diffusion equation (13) can still be used to obtain 
the average power transferred by the wave to the resonant electrons. The evolution of 
electron kinetic energy is given by 


dE c 

dt 


d 1 2 f df 1 2 

— {-n e m e v ) = / av— — n e m e v . 

dV 2 1 J dt 2 


(14) 


Using (13) and one integration by part yields 

dEc _ nn e e 2 \E k)U \ 2 uj df 
dt m e k 2 dv 

consistent with electron acceleration when df /dv < 0 or deceleration when df /dv > 0 . 

Moreover, conservation of the energy given by d{\E k ^\ 2 /^ + E c )/dt = 0, yields the 

following equation, d\E k ^\ 2 /dt = 2'y\E k)U \ 2 , describing time evolution of a given Fourier 

component of the electric held fluctuations. When these electric held fluctuations are 

produced by Langmuir waves with u ~ oj pe , the coefficient 7 corresponding to Landau 

damping or growth of the waves is 



7 = 


df_ 
2 k 2 dv 


V=OJpe/k 


(16) 


Summing over all Fourier components of the wave held, the evolution of the spectral 
energy density is also given by 


dW(k,t) 

dt 

March 6, 


= 2 7 LF(M)- 

2015, 2:35am 


(17) 
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Equations (12) and (17), are the standard quasilinear diffusion equations describing the 
self-consistent evolution of a population of electrons and a spectrum of Langmuir waves 
these electrons may generate [Pines and Schrieffer, 1962; Vedenov and Velikhov, 1963; 
Drummond and Pines, 1964], The coefficient 7 , giving the rate of growth or damping 
of each spectral component of the Langmuir wave electric field, is an average rate of 
stimulated emission and absorption of the waves by the particles, a point which can 
be made more transparent in the semi-classical derivation of these quasilinear diffusion 
equations [Pines and Schrieffer, 1962; Melrose, 1980; Tsytovich, 1995]. The quasilinear 
diffusion equations generally include terms describing the spontaneous emission of waves 
and the dynamic friction on the particles. 


3. Broadening of the wave-particle resonances by particle and wave scattering 

The basic rates of wave absorption or emission corresponding to particle acceleration 
or deceleration are affected by the presence of background scattering of particles or waves 
in the medium. In the following, the effect of scattering is taken into account in the 
quasilinear diffusion equations as a broadening of the wave-particle resonance. Indeed, 
the resonant nature of the interactions can be affected by any physical mechanism which is 
able to destroy coherence, say on a time scale Ac v, resulting in broadening of the resonance 
function, for instance into a Lorentzian form: 


tt5(uj — kv) —> 


Ac v 


(u — kv) 2 + Acu 2 


(18) 


Such Lorentzian form can be obtained by keeping v finite in Equation (10), making it 
clear the physical interpretation of resonance broadening due to wave damping, 


Acu = v, 


(19) 
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i.e., the resonant interaction between an electron and a Langmuir wave is limited by the 

life-time of the wave. Formally, resonance broadening consists in the substitution 

7t<5(u; — kv) —> R(lo — kv, Au), (20) 

with R not necessarily a Lorentzian, in the quasilinear diffusion equations. It is here 
interesting to recall a point made by Van Vleck and Margenau [1949] in a related con¬ 
text, which is that line-broadening bridges the gap between resonant and non-resonant 
absorption processes as the level of scattering in the medium increases. Such distinction 
between resonant and non-resonant acceleration by waves, involving resonance-broadening 
as a transition, is also at the core of a recent classification scheme recently proposed by 
Bian et al. [2012] for stochastic acceleration during flares. 

In the following section, we discuss two important sources of resonance broadening in 
the beam-plasma system. These are the scattering of beam electrons and the scattering 
of beam-generated Langmuir waves in turbulent plasmas. 


3.1. Resonance broadening due to particle scattering 

One main effect of background scattering of the particles on the wave-particle interac¬ 
tions can be understood by modelling such scattering as a additional random perturbation 
in the original equations of motion : 


dv e dx w . . . 

H = m e E(X ' t) ’ M =V + m ’ (21) 

where ((t) is a Gaussian white noise. Scattering enters here as an external source of 
stochasticity in the equations of motion for the particles. Notice that such scattering 
does not produce any change in the kinetic energy of the particles, acting also in ab¬ 


sence of the electric field. We shall see, however, that scattering affects the particle 
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acceleration/deceleration rate due to the electric field and, from energy conservation, also 
modifies the rate of damping/growth of the electric energy. In fact, scattering results in 
broadening of the wave-particle resonances because it interrupts at random the coherent 
emission/absorption of the waves by the particles. In this case, one has to deal with the 
random nature of the integral / 0 °° dt(e lkx ~ lwt ), defining the resonance function and, hence, 
the particle path is written in the form 

x = vt + Ax, (22) 

where the perturbation Ax accounts for the random scattering of the particles around 
their free-streaming trajectory. Under the additional assumption that Ax is a Gaussian 
random process with (Ax 2 ) oc t a , the resonance function becomes broadened around 

c o — kv , in the form of a Lorentzian for a = 1, of a Gaussian function for a = 2 or in the 

form of a Airy function for a = 3. 

3.1.1. Spatial diffusion 
The case a = 1, with 

(Ax 2 ) = 2 Kt, (23) 

represents a standard spatial diffusion of the particles around the free-streaming trajectory, 
with a spatial diffusion coefficient n. For a Gaussian process £ with zero average, (£) = 0, 
and finite variance (£ 2 ) / Oa cumulant expansion [Kubo, 1962] provides the relation 

( e «) = e <« 2 > /2 , (24) 

where the brackets denote averaging over the probability distribution function (PDF) of 
£, here assumed to be Gaussian. There is a-priori no need to restrict to Gaussian PDFs 
except as a simplifying assumption. Taking £ = ikAx, we obtain that the resonance 
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function is given by 

fOO A / j 

R(u - kv, Au) = / dte ikvt ~ iut e~ Aujt = --—--— , (25) 

v ’ Jo (u - kv) 2 + Au 2 y ’ 

where the resonance width Au is given by the diffusive time-scale over one wavelength of 

the electric field, corresponding to 

Au = nk 2 . (26) 


3.1.2. Velocity-space diffusion 

Equivalently, the perturbations in the particle path can be written in terms of velocity 
fluctuations 

x — (v + A v)t, (27) 

and, furthermore, we assume that the fluctuations Av evolve diffusively, 

(Au 2 ) = 2 Dt, (28) 

where D is the diffusion coefficient in velocity space. As a consequence, the resonance 
function becomes broadened around c o — kv as 

pOO o 7T 

R(u - kv, Au) = / dte ikvt - iujt e~ (Au}t) = — C\i{u - kv)/Au\. (29) 

Jo Au 

This expression involves a Airy function C(z) which satisfies the equation d 2 C/dz 2 — zC = 
7T _ 1 [Gradshteyn and Ryzhik , 1980]. The resonance width is given by 

Au = ( k 2 D) l/ 3 . (30) 

Let us provide few examples of possible sources of velocity-space diffusion. In the original 
work by Dupree [1966], the source of velocity-space diffusion is the wave electric field itself. 
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This means that the resonance function (29) involved in the expression for the velocity 
space diffusion D given by Eq.(8) is now a function of D, hence we can rewrite (8) in the 
form 

D = ^-\ [ dkduS E (k,u) f°° dte ikvt ~ iu}t e~ k2Dt3 . (31) 

rnj J J Jo 

This is a non-linear integral equation for D. If we take the limit of weak fluctuations 
amplitude, D —> 0 in evaluating the resonance function R then (31) reduces to the familiar 
quasilinear result (11) with a sharp resonance function , i.e. R{u — kv,Au —* 0) —> 
7tS(co — kv). In other words, the quasilinear equations are modified by broadening as 
the strength of the electric held fluctuations becomes large. Physically, the resonance 
broadening effect corresponds to secular diffusion of particle trajectories which results 
from the interaction of the particles with the turbulent electric held fluctuations. We 
notice that a Gaussian velocity space diffusion with (Av 2 ) oc t leads to Gaussian spatial 
super-diffusion with (Ax 2 ) oc f 3 , hence the form of the resonance function (29). 

The source of velocity-space diffusion involved in the broadening of the wave-particle 
resonance may also be unrelated to the wave electric held but instead due to background 
scattering of the particles for instance at constant kinetic energy. Beam electrons strongly 
interact with the ambient magnetic held as the latter imposes a direction to their propaga¬ 
tion. Moreover, magnetic fluctuations around the mean ambient held produce pitch-angle 
scattering which acts as an additional source of velocity space diffusion and, hence, of res¬ 
onance broadening. Pitch-angle scattering results in fluctuations in the partition between 
translational (along the ambient held) and rotational energy of the particles in a mag¬ 
netized plasma and when this is due to a spectrum of low-frequency magnetic perturba¬ 
tions, pitch-angle scattering occurs at constant kinetic energy E c . A rigourous treatment 
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of pitch-angle scattering would require a 3D analysis, however let us emphasize few main 
aspects of it in a ID model. Let us still call v the velocity parallel to the magnetic held, 
V = \J2E c /m e the electron speed and // = v/V the pitch-angle cosine. Interactions of the 
electrons with magnetic held fluctuations produce increments of pitch-angle cosine A fi at 
constant V, and hence increments of parallel velocity Av = VA/j. Squaring the result 
and averaging leads to 

(A,, 2 ) = ^(A,, 2 ). (32) 

m e 

Moreover, let us assume that pitch-angle scattering of the particles produces pitch-angle 
diffusion with 

(A/j, 2 ) = 2 D^t, (33) 


where is the pitch-angle diffusion coefficient. Hence, pitch angle scattering acts as a 
source of diffusion in velocity space but at constant kinetic energy E c , i.e. 


(An 2 ) = 2 


2E c DM 

—-- 1 

m e 


(34) 


a relation which is analogous to (28). And therefore, the resonance function becomes 
broadened with a width given by (30) where D has now the meaning of the coefficient of 


diffusion for the parallel velocity at constant energy, which according to (34) tantamount 
to replacing D —> (2 E c /m e )D^ in (30), leading to 



which is the resonance width associated with such scattering mechanism. The above 
reasoning is clue to Volk [1973] in his treatment of resonance broadening by pitch-angle 
scattering. Resonance broadening due to collisions, D ijli taking its collisional value value 
in (35) is discussed by Stix [1992], 
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3.2. Resonance broadening due to wave scattering 

While resonance broadening due to scattering of particles has been widely discussed 
in the past, there is also a possible source of resonance broadening, which is not due to 
particle scattering but which is instead due to wave scattering. 

3.2.1. Wave-number diffusion 

For clarity, we first discuss the effect of wave scattering through the following model, 
which is a variant of the Kubo-Anderson oscillator [ Anderson , 1954; Kubo, 1954], The 
model describes the motion of electrons, and their stochastic acceleration or deceleration, 
in a plane wave electric field whose wave-number is a random function of time: 


dv 

dt 


— (4, u 

77 ) _ \ 


0 i[(k+Ak)x—cjt] 


+ CC , 



where, for simplicity, ((t) is taken to be a Gaussian white noise. According to (8) the 
expression for the diffusion coefficient in velocity space is given by 

p2 ^ r-OO 

D = — \E k J 2 dte ikvt - iu}t (e iAkvt ) . (37) 

mj ’ Jo 

The evolution of the wave-number in this model is a standard diffusion : 


(A k 2 ) = 2 D*t, (38) 

with D* the diffusion coefficient in wave-number space. Therefore, we obtain that the 
diffusion coefficient in velocity space D is related to the diffusion coefficient in wave- 
number space D* by 

D = — \E ku \ 2 / dte ikvt ~ iwt e~ D * v 1 . (39) 

mj ’ Jo 

We have implicitly assumed that the effect of particles scattering is negligibly small com¬ 
pared to that of the waves. Equivalently, the diffusion coefficient in velocity space can be 
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written in the form 

e 2 ~ 

D =—\E k u \ 2 R(uj - kv,Au) , ( 40 ) 

ml 

with the resonance function still involving the Airy function given by (29) but the reso¬ 
nance width is now given by 

Acu = (L>V) 1/3 , ( 41 ) 

with D* the diffusion coefficient in k- space. While the plane wave is randomly scattered 
it also undergoes Landau damping when accelerating the particles. The average power 
transferred to the particles by a wave undergoing random change in k leads to the coef¬ 
ficient of Landau damping of this wave which can be calculated in a manner similar to 
what was done for obtaining the Landau damping rate (16) of an unscattered wave. The 
result is 

7 = ~ kv i Aw), ( 42 ) 

which reduces to the standard expression for Landau damping of a Langmuir wave given 
by (16) when u ~ c o pe and when the level of wave scattering becomes small, i.e. D* —> 0 
or Acn — > 0. 

Wave-number diffusion has been widely discussed in the context of the modelling of the 
propagation of Langmuir waves in a plasma with background density fluctuations. Above, 
we have shown that such wave-number diffusion leads to broadening of the wave-particle 
resonances. In the following, we derive a system of resonance-broadened quasilinear dif¬ 
fusion equations which describes the relaxation of an electron beam in a plasma with 
ambient density fluctuations. 
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3.2.2. Broadening of the wave-particle resonances in a plasma with density 
fluctuations 

The dynamics of both electrons and Langmuir waves are still taken along x, the direction 
of the external magnetic field. The plasma density is written as n e [l + n(x, £)] with n e the 
constant background density and n(x, t) the relative density fluctuation, which is assumed 
to be weak, i.e. n(x,t) <C 1. In general, the characteristic wavenumber q of low-frequency 
density fluctuations is much smaller than the characteristic wavenumber k associated 
with high-frequency Langmuir waves, so we can make the WKB approximation and treat 
the Langmuir waves as plasmons or quasi-particles. Ambient density fluctuations are 
associated with a change in the local refractive index experienced by the waves, so that 
wave scattering can be modeled by the equations of motion for the Langmuir plasmons, 
which are given by 

dk 1 dn 

Tt = -r^^ F{x ' t) ' (43) 

dx 
dt 

where F(x,t) is the ’’refraction force” acting on the quasi-particles, u> pe = \J 47rn e e 2 /m e 
is the local plasma frequency and v g = 3 v^ e k/uj pe is the group velocity of the quasi¬ 
particles. From Equation (43), we see that random refraction induces stochastic change in 
wavenumber, resulting in a diffusion of the wave energy density in A;-space. An expression 
for the diffusion coefficient of wave energy in /e-space is thus derived, in the very same 
quasilinear approximation, 

2 2 

uj 7T r°° r 00 

D* = dq dn q 2 S n (q , - qv g ), (45) 

Z . 7—00 . 7—00 
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where by definition (n 2 ) = f^dqf^dU S n (q,kl). In the particular case where density 
fluctuations are due to low-frequency compressive waves with a specific dispersion relation 
n = Q(q), then S n (q, ft) = S n (q)S(Q - f1(g)) and D* = (w 2 e 7T 2 /2) f™ 00 dqq 2 S n (q)5(Sl(q) - 
qv g ). A set of quasilinear equations describing the interaction between Langmuir waves 
and beam electrons in the presence of background density fluctuations is therefore, 


dt 

dW(k,t) _ ttw 3 


4n 2 e 2 d W(k = 

: Upe/V,t) df(v,t ) 

(46) 

ml dv 

v dv 


d dW(k,t ) 

„ dk dk 

(47) 


dt k 2 

which includes a wave-number diffusion in the kinetic equation for the plasmons. As a 
consequence, this system is equivalent to the following quasilinear diffusion equations with 
a broadened resonance function : 


df(v,t) 4yre 2 d f^^ rn , A , df(v,t) 

~dT = ~^Tv J dkW(k ' t)R{ - kv ’ - 

dW(k,t ) uj 3 f df(v, t) 

- Ql - = -gi J dvk—^—R{u pe - kv, Aui)W(k, t), 

where the resonance function is given by (29) and the resonance width is 

Aw = (L>V) 1/3 , 

with D* given by (45). 


(48) 

(49) 


(50) 


4. Generalization to three dimensions 

Let us consider the generalization of the quasilinear diffusion equations to three dimen¬ 
sions including a broadened resonance function. These are 
df(v,t) Ane 2 d 


d k ^lV k fl ( ^- k . V ,A W ) S/(v '*) 


dt ml di>i J k 2 

describing the diffusion of resonant electrons in velocity space, and 

k df(v,t 


dvj 
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~W = UJ ^ 


W&rjfo . 2(A^,’ fShpi 


(51) 
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giving the rate of growth or damping of spectral energy in the Langmuir waves. They 

reduce to the standard 3D quasilinear diffusion equations when the broadening width goes 
to zero, i.e. R{u pe — k.v, Au —* 0) —* 7rS(u pe — k.v). The resonance function appearing 
in these equations is defined as 

R{u pe - k.v, Au) = f dt{e ik - vt ~ iuipet ). (53) 

Let us first consider the effect of random perturbations Av in the electron velocity, i.e. 
v —)■ v + Av, hence 


R(u pe — k.v, Aca) = / dte 


ik..vt—iuj ve t 


(e 


ik.Avt^ 


= / dte 


ik.vt— iujpet^— ((k.Av) 2 )t 2 /2 


(54) 


Introducing the diffusion tensor in velocity space, which is defined via 


(AviAvj) = 2Dijt, (55) 

we obtain from (54) the following expression for the resonance function 

R(co pe - k.v, Au) = j dt^-^e-^ 3 , (56) 

where the resonance width Au is given by 


Au = Dijkikj. 


(57) 


Summation over repeated indices is implicitly assumed in the previous expression. In the 
same way, let us consider random perturbations Ak in the wave-vector, i.e. k —)■ k + Ak, 
hence 


R(u pe — k.v, Au) = / dte 


ik..vt—iu) pe t 


(e 


zAk.vt^ 


= / dte 


zk.vt— iujpet^— ((Ak.v) 2 )r/2 


(58) 


Introducing the diffusion tensor in wave-vector space, 


(59) 
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we can again write the resonance function in the form of a Airy function given by Eq.(56) 

but with the resonance width Acn given by 


Ac jj = D^ViVj. 


(60) 


Broadening of the wave-particle resonance occurs here as a result of the random refraction 
of beam-generated Langmuir waves as they propagate in fluctuating three-dimensional 
density fluctuations, and therefore, 


dk 

dt 



dx 

dt 


= v 


9’ 


(61) 

(62) 


where n = n/n 0 is the relative level of density fluctuations and v s = 3\ 2 De u> pe k is the group 
velocity. The resonance width (60) involves the diffusion tensor D* ;) in wave-number space, 
which is given by 


2 

D ij = j dq J q,fi)5(n - q.v ff ), (63) 

where by definition (h 2 ) = f^dqf^dfl SV^q, hi), meaning that A n (q, hi) is the spec¬ 
trum of (relative) density fluctuations. We here assume a given level of large-scale density 
fluctuations of an unspecified origin. The density fluctuations do not have to be pro¬ 
duced by waves in which case hi and q are not related by a specific dispersion relation 
Q(q). However, when the background density fluctuations are produced by low-frequency 
compressive waves, self-consistency dictates that energy change in the Langmuir quasipar¬ 
ticles is accompanied by energy change in the low-frequency compressive waves resulting 
in Landau damping/growth of the low-frequency modes on the high-frequency modes, 
the latter behaving as quasi-particles [Vedenov et al., 1967]. In order to respect energy 
conservation, the system must be closed by a third equation for the level of compres- 
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sive density fluctuations which results in D*j becoming a dynamical variable, and, hence, 
also the resonance width Aca, with their time-dependence adjusting self-consistently to 
the non-linear channeling of the energy between the various modes of motion. In other 
words, the resonance broadening effect due to Langmuir wave scattering by low-frequency 
compressive waves is essentially a manifestation of mode coupling. In the following we 
focus on mode-coupling to ion sound waves but Langmuir wave scattering may be due 
to other compressive modes, such as kinetic Alfven waves [Bian et a/., 2010; Bian and 
Kontar, 2010]. The spectral properties of these electromagnetic compressive modes are 
well-documented in the solar wind [Bale et al, 2005; Mozer and Chen , 2013]. 

4.1. Angular scattering only 

Before continuing to the topic of resonance broadening by mode-coupling, let us discuss 
briefly the problem involved, in general, in determining the shape of the resonance function 
and, in particular, in the case where there is only angular scattering of the particles or only 
angular scattering of the waves. The resonance function entering the quasilinear diffusion 
equations involves the quantity c o pe — k.v = co pe — kv [sin 9 sin d'(cos <f> cos qb' + sin (f) sin 0') + 
cos 9 cos 9'], which is written here in the spherical coordinate systems (v, 9, <fi) and (k, 9', <p') 
in velocity space and wave-number space, respectively. An essential problem arises when 
trying to determine the shape of the resonance function when diffusion occurs only in 
angle, 9 for the particles or 9' for the waves. The problem is that none of these diffusion 
processes are characterized by Gaussian PDFs. This is quite obvious by noticing that 
these PDfs must be periodic functions of the angles involved and, therefore, cannot be 
Gaussian. There are four angles involved in the problem and they can be reduced to two 
by restricting our considerations to not too large angular spread of the particles or the 
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waves around 6 ~ 0 or O' ~ 0 corresponding to the direction of the mean held, i.e. 

g 2 012 

uj pe — k.v = c o pe — kvcosOcosO' = L0 pe — kv( 1 ——)(1 - —). (64) 


Let us first consider angular scattering of the particles assuming a narrow distribution of 
waves having O' = 0, such that 

02 

u pe - k.v = u pe - kv(l - —). (65) 


Pitch-angle diffusion of the particles results in that P(8,t), the PDF for the pitch-angle 
6 evolves according to the diffusion equation 


dmt ) _19 dP(0,t) 

dt odo 6Dee do 


( 66 ) 


when sin 6 ~ 0. Here, Dgg is the component of the diffusion tensor corresponding to 
particle scattering in 6 only. The solution of this equation can be approximated by a 
Gaussian only for 6 ~ 0 and when Dggt <C 1, in which case 


< 6 2 >= 2D w t. 


(67) 


This the domain of validity of the approximate approach discussed in Section III above, 
for the treatment in ID of resonance broadening by pitch-angle scattering of the particles, 
where we had introduced instead of 6 the notation /j = cos 6 for the pitch-angle cosine of 
the particles. Outside of this domain of validity, analytical determination of the resonance 
function must rely on the more accurate PDF P(/i, t) which is known to be given by 
series expansion in terms of Legendre polynomials. The averaging over P(n,t) which 
is involved in the determination of the resonance function becomes a more complicated 
analytical procedure, probably amenable to numerical solution only. However, the above 
approximate treatment provided us with an estimate for the resonance width, despite that 


the precise ’’line-shape” may not well be the correct one. 
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The perfectly symmetric situation of angular scattering of the waves at constant mod¬ 
ulus k requires similar considerations, as in this case 

0i'2 

U pe - k.v = UJ pe — kv(l — —), (68) 

and angular diffusion of the waves results in that P(0', t), the PDF for the angle O' evolves 
according to the diffusion relation 

< 0 2 >= 2 D* g , e ,t, (69) 

where Dq, 0 , is the component of the diffusion tensor corresponding to wave scattering in 
O' only. Treatment of the simultaneous scattering of waves and particles is also not quite 
a simple task as this requires determination of the joint PDF P(0,0',t) in order to carry 
the averaging procedure with respects to these random variables. 


5. Resonance broadening due to mode-coupling of Lamgmuir waves with ion 

sound waves and weak-turbulence theory 

Let us focus on the specific situation of coupling to the ion-sound waves described by 

the Zakharov [1972] equations 


j_dE_ 

^pe (P 


AL V 2 E = 


n 0 mi 


' d ' 2 

__ r 2 v 2 

dt 2 s 


n 


n 0 


n 

2 % 
1 


E, 




47T 


(70) 

(71) 


where c s is the sound speed and \r> e is the Debye length. The first equation is a 
Schrodinger equation describing the evolution of the Langmuir wave electric field in low- 
frequency density inhomogeneities. The second equation describes the dynamical evo¬ 
lution of these density inhomogeneities when the latter are assumed to be produced by 
sound waves. These wave equations are coupled through the terms on their right-hand 
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side. These terms describe refraction and ponderomotive effects, e.g. [ Thornhill and 
Ter Haar, 1978; Robinson, 1997]. The Zakharov equations exploit the separation of time 
scales between the dynamics of Langmuir and ion-sound waves, so that the ponderomotive 
force represents in fact the time average back-reaction of the high-frequency oscillations 
on the low-frequency density inhomogeneities. The ponderomotive force is the spatial 
gradient of a potential, i.e. f = VI/ with U = |i?| 2 /87mo, and this force can be shown 
to act primarily on the electrons. In addition to the separation of time-scales inherent 
to the Zakaharov equations, one can further exploit a separation of length-scales by con¬ 
sidering the particular case of non-local three-wave interactions involving an ion-sound 
mode having a wave-vector and frequency given by [q, f2(q)] and two Langmuir modes, 
[k, <u(k)] and [k + 5k, cu(k + 5k)], under the restriction 5k <C k for non-local mode cou¬ 
pling. Therefore, the two resonance conditions for mode coupling are q + k = k + 5k 
and 11 (q) + o;(k) = cu(k + 5k). By Taylor expanding the second condition, we obtain 
that 12(q) = 5k.VA,o;(k) and using the first one, which is q = 5k, the resonance condi¬ 
tions for the non-local three-wave interactions are found to reduce to only one condition 
ff(q) = q.v 3 , i.e. ±qc s = 3A|, e a;p e q.k, which is similar to a wave-particle resonance 
condition cn(k) = k.v. Non-local resonant interactions between the large-scale and the 
small-scale waves (quasi-particles) are therefore expected to mediate local diffusive energy 
transfer among the small-scale waves in analogy with wave-particle interactions. By local 
diffusive energy transfer we mean that the spread of spectral energy associated with the 
small-scale waves involves small steps in wave-number space, having 5k <C k and that this 
spread is a diffusion. Since the first Zakharov equation has the form of a Schrodinger equa¬ 
tion describing the dynamics of Langmuir waves in the a slowly varying potential n/no, 
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it is standard to apply WKB analysis to it. In order to do that it is convenient to use 
a generalization of the Fourier transform to spatially inhomogenous systems [Gershgorin 
et al, 2009]. The window transform of E(x,t) is defined by 

r[£(x)] = E(x, k) = J dx 0 w(e*|x - xo|)^(x 0 )e* k - Xo , (72) 


which is a kind of wavelet transform that in the case where the window function 
w(x) = exp(— x 2 ) is called a Gabor transform. The parameter 1/e* is the width of the 
window which is taken to be much smaller than the characteristic length of the inhomo¬ 
geneity but much larger than the wavelength of the Langmuir waves that propagate in 
the inhomogeneous medium. Note that the Gabor transform can be viewed as a localized 
Fourier transform around x with support 1/e* : when e* goes to zero, the filtering width 
goes to infinity and the Gabor transform becomes a Fourier transform. Applying this 


transform to Equation (70), we have 


9E kx 

dt 


- * £ 


Ak.x + *V I —— 

\2n 0 


<9E k , x 
‘ <9k 


3AnJk.Vi9 k . ; 


Multiplying by E k x , the imaginary part of (73) gives the Liouvillc equation for the small- 
scale wave energy density. Information on the phases of the small-scale waves is now lost 
as these waves are treated as quasi-particles. Moreover, let us introduce the wave-action 
density, defined as 


(k)’ 


and the group velocity w g = 3A|, e a;p e k, therefore, the WKB limit of the Zakharov equations 
are written as in Vedenov et al. [1967]: 


d N k 
dt 


+ v w 


d N k 
<9x 



V 



dN k 
' dk 


0, 


(75) 
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nom ,g-e^](A) = _i_ v2 / (jk ^ (76) 

The Zakharov [1972] equation were first written by Vedenov et al. [1967] in this WKB 
form, which corresponds to restricting the consideration to non-local three-wave coupling 
with kL,kL' ~ k 3> ks = q in the process L + S = L'. This system is similar to the 
field/particle system of Eqs.(l)-(2) we started from, so that quasilinear diffusion equations 
can be derived from it. Linearizing the first equation (75) with respect to AT = N^ + SN^ 
i(fi(q) — q.v g )5N k = (uj pe n/2n 0 ) q.dN^/dk, and substituting 5N k in the second equation, 
one obtains the dispersion relation connecting the frequency fl(q) and the wave-vector q 
of the ion-acoustic wave which now includes an imaginary part 7 q . By considering the 
imaginary part y q of f2(q) to be small compared with the frequency fi(q), i.e. y q <C fl(q) 
then 


™p e q 
8 n 0 mi.c s 


fi(q) 


where in the last expression fi(q) = ±qc cs is used as dispersion relation. The imaginary 
part 7 q is the coefficient describing the Landau damping or growth of the ion-acoustic 
waves on the Langmuir plasmons. When it is positive the coefficient (77) is the growth 
rate of the so-called modulational instability [ Thornhill and Ter Haar, 1978]. 

Moreover, we have seen above that Langmuir plasmons undergo a diffusion in wavenum¬ 
ber space in the presence of a spectrum of density fluctuations. Here, this diffusion pro¬ 
cess in wave-number space occurs as a result of the action of the stochastic refraction 
force of the sound modes on the qnasi-particles. The diffusion coefficient is given by 
D*j = (T0J pe /4) f dqq i qjS n {q)5(Tl(q) — q.v 9 ). Therefore, in the non-local limit and under 


DRAFT 


March 6, 2015, 2:35am 


DRAFT 



BIAN ET AL.: RESONANCE BROADENING IN QUASI-LINEAR RELAXATION X - 29 

the random phase approximation, the Zakharov equations take the common quasilinear 
diffusive form given by 


dN{ k) d dN{ k) 
dt dh ij dkj 


dS n ( q) 
dt 


7q'S'n(q)- 


(78) 

(79) 


On the other hand, we have seen that wave-number diffusion is a source of broadening 
of the wave-particle resonance. Therefore, the quasilinear diffusion equations describing 
the interactions between Langmuir waves and electrons will involves a resonance func¬ 
tion R(u — k.v, Acn) which is broadened due to the mode-coupling of Langmuir waves 
with sound waves described by Eqs.(78)-(79). As shown in Section 4, the wave-particle 
resonance width is given by Acn = {D*-ViVj) 1 ^. 

To summarize, we have obtained the following resonance-broadened model of weak- 
turbulcnce: 


df(y) = 
dt 

dN( k) 
dt 

dS n ( q) 
dt 

where 


^ n df(v) 
dvi lj dvj 

(80) 

= 7kiV(k), 

(81) 

= 7q^n(q), 

(82) 


n 4vr e 2 uj pe 

JJij — 


kjkj 


mi 


7k = / dv 


N(k)R(ujp e — k.v, Acn), 
Riojpe — k.v, Acn), 


k 2 

k df{v) 


7q = 


™ 2 ve q 

8 n 0 miC s 


k 2 dv 

I _q v ») 


(83) 

(84) 

(85) 


and where the wave-particle resonance width entering R(u pe — k.v, Acn) is given by 


Acn = (D^ViVj) 1 ^, 
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which is related to the dynamical level of sound waves fluctuations through 

_ 2 

D *j = J dqg i g j 5 n (q)5(Q(q) - q.vj, (87) 

where T2(q) = ±gc s . In this formulation of the Vedenov et al. [1967] equations, the effect of 
mode-coupling between the Langmuir plasmons and the sound-waves is taken into account 
through a broadening of the wave-particle resonances. 

Notice that the level of sound waves *57 (q) entering (87), which is also involved in 
the expression for the broadening width (86), is depending on time via Eq.(82) which 
involves y q , and hence dN^/d k . This dynamical broadening of the wave-particle in¬ 
teraction reflects the self-consistent level of Langmuir wave scattering by an evolving 
spectrum of sound modes in the beam-plasma system. Disregarding time evolution of 
the sound waves given by Eq.(82) results in Eqs.(80)-(81), together with Eqs.(83)-(84) 
and with Eqs.(86)-(87), being a description of the relaxation of an electron beam in a 
prescribed level of density fluctuations unaffected by the Langmuir waves. These den¬ 
sity fluctuations may not necessarily be produced by sound waves for an arbitrary dis¬ 
persion relation T2(q). These density fluctuations may not necessarily be produced by 
waves when no such dispersion relation exists, in which case we formally need to take 
D*j = 7ru^ e /4fdq J'dnqppSJq, Q)S(Q - q.v g ) instead of (87) in the expression for the 
wave-particle resonance width (86). 

6. Summary and discussion 

In general, the resonance function of quasilinear theory has a finite width and can be 
written in the form given by 

r oo 

R(u-kv, Au)= / dt(e ikx - iut ), (88) 

Jo 
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where Aca is the resonance width. The quasilinear diffusion coefficient then takes the form 

D = —— dk du Ssik, ca) R(uj — kv, Aw), 
mj J J 

where Sfg^/qca) is the spectrum of electric held fluctuations. We are interested in electric 
held fluctuations which are produced by Langmuir waves whose dispersion relation is 
u>(k) ~ uj pe = \J 47rn e e 2 / m e , and hence we can write Ssik^u ) = S^/c) <5(ca — ca pe ). The 
spectral energy density of Langmuir waves is given by W{k ) = 2(SE(k)/8n ), where the 
factor of two accounts for the electric energy and the kinetic energy of the thermal electrons 
participating in the wave motion. Hence, the quasilinear diffusion coefficient is given by 

D-— —— J dkW(k)R(uj pe — kv, Alu) . (90) 

By conservation laws, statistical acceleration of the particles by the waves, or their decel¬ 
eration, means damping or growth of the waves. Therefore, the rate of Landau damping 
or growth of Langmuir waves is given by 

7 = J dvk-^-R(uj pe - kv, Au). (91) 

The quasilinear diffusion equations, involving a broad resonance function, have the stan¬ 
dard form given by the system 



df(v) = d df(v) 
dt dv dv 


dW(k) 
dt 


2'yW{k), 


(92) 

(93) 


which by construction respects conservation of the energy. In practice, it is useful to 


approximate the broad resonance function, entering the quasilinear diffusion equations, 
by a square function 


R[ojpe — kv, Aca) 


7T 

Alu ’ 


o, 


uj pe — kv | < Aca, 
t o pe — kv | > Aca, 
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having the property that i?(w pe — kv , Aw) n S(u pe — kv ), the sharp resonance being 

obtained from x = vt and w = co pe in (88) . 

Turbulent scattering of particles results in velocity fluctuations Av around the free- 
streaming trajectory so that we have instead x — (v + Av)t in (88). When these fluctua¬ 
tions evolve diffusively, (An 2 ) = 2 Dt, then the width of the resonance is given by 

Aw = (Dk 2 ) 1/3 . (95) 

In Dupree [1966] work D is the quasilinear diffusion coefficient itself, so that Aw is a func¬ 
tion of the electric field amplitude. Pitch-angle scattering off turbulent magnetic field fluc¬ 
tuations can also broaden the wave-particle resonance leading to Aw = (2E C D^k 2 /m e )i/ 3 
which tantamount to replacing D —y (2 E c /m e )D^> in (95), with D T) the pitch-angle dif¬ 
fusion coefficient. Resonance broadening due to collisional pitch-angle scattering, 
taking its collisional value, is discussed by Stix [1992], Langmuir wave scattering off back¬ 
ground density fluctuations can also result in broadening of the resonance. In particular, 
when the beam-driven Langmuir wave undergoes stochastic refraction, its wave-number 
evolves diffusively so that k —> k + Ak in (88). Assuming diffusive evolution in wave- 
number space (Ak 2 ) = 2 D*t, the resonance width is then given by 

Aw = (L>V) 1/3 , (96) 


where D* is the diffusion coefficient in wave-number space. Finally, let us discuss the 
consequence of resonance broadening on the development of the beam instability. From 
the broadening frequency Aw, one can also define a scale Av in velocity-space given by 


Av = 


Aw 

IP 


(97) 
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Using the square as a representation of the resonance function, the growth rate of the 

beam-plasma instability can also be written in terms of Av, as 


j(Av) ~ 


nu pe fjupe/k + An/2) - f(upe/k - An/2) 

2 k 2 Av 


™pe A f 
2 k 2 Av 


(98) 


In other words, the basic time-scale 7 for the development of the beam instability is 
now a function of the resonance width Au through Av. Resonance broadening filters 
out variations in the electron phase-space distribution which are smaller than Av. As 
phase-space derivatives also yield the rate of growth of the beam instability in Eq.(98), 
7 will in turn depend on which scale the phase-space has been filtered by the underlying 
broadening processes. Notice that when Av —> 0, the expression for the growth rate 
Eq.(98) including broadening reduces to the standard Landau damping/growth rate, i.e. 
7 (Au -> 0 ) -► (irujp e /2k 2 )df /dv\ v=UJpe /k- This means that for a given distribution function 
f{y) including a beam of extent A Vb, the growth rate y(Av) can either be positive or 
negative depending on the value of the filtering scale Av in velocity space as illustrated 
below. Let us consider qualitatively the situation of the beam instability depicted in 
Fig 1. Due to the presence of the electron beam, Langmuir waves are expected to grow 
with a phase speed of the order of the beam speed ~ 17 . However, the presence of 
background scattering has the effect on broadening the resonant region in velocity space. 
For small resonance width Av = V\ — v 2 <C Au&, the average slope of electron phase-space 
distribution is still positive, therefore the particles loose energy and the waves grow but 
at a reduced rate. The effective growth rate is given by Equation (98) for 'y(Av). For 
a large resonance width i.e. when Av = v ?J — iq 3> A Vb, the average slope becomes 
positive and the waves damp as now -/(Av) < 0. In the situation of Fig (1), the condition 
Av A Vb sets the regime of the suppression of the beam instability due to the presence of 
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background scattering of the particles or of the waves. The growth rate (98) turns out to 
be negative despite the existence of a region of positive slope for the distribution function, 
because small-scale features in velocity space, smaller than An, and which should act as a 
reservoir of free-energy for the Langmuir waves to grow, have been effectively filtered out. 
Resonance broadening is a general mechanism at the origin of the weakening and possible 
suppression of the beam-plasma instability. As a consequence, it can be used as a tool 
to explain stable inverted populations of fast electrons commonly observed in-situ in the 
solar-wind (P. Kellogg, private communication). We refer to the work by Lin et al. [1986], 
where electron distribution functions are observed not to exhibit the strong plateauing 
predicted by quasi-linear models. 

With respect to Langmuir waves scattering off background density fluctuations, we 
obtain the following criterium for the suppression of the beam instability : 


D *>-f 


uL ( Av h x 3 


Vb \ V b 


(99) 


When needed, a more accurate criterium can be obtained by directly solving the marginal 
stability 7 (Acn) = 0 condition, Equation (98). 

To summarize, we have shown that resonance broadening provide a unified framework 
accounting for both the effects of particle and wave scattering during the qnasilinear 
relaxation of electron beams in turbulent plasmas. Resonance broadening is essentially 
equivalent to filtering-out small scale features of the electron phase-space distribution. 
As a beam is such a small scale feature in velocity space, resonance broadening provides 
an intuitive picture for the stabilization of the beam instability. Starting directly from 
the electrostatic Zakharov equations, we derived a resonance broadened version of weak- 
turbulence theory, where broadening of the wave-particle resonance occurs as a result of 


DRAFT 


March 6, 2015, 2:35am 


DRAFT 



BIAN ET AL.: RESONANCE BROADENING IN QUASI-LINEAR RELAXATION 


X - 35 


non-local mode coupling to ion-sound waves. In this case, broadening is non-linear and 
dynamical reflecting the self-consistent level of ion-sound waves generated by Langmuir 
waves in the system. In the future, we plan to extend this formalism to account for the 
mode-coupling of beam generated Langmuir waves to kinetic Alfven waves. 
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Figure 1. An example of weakening and possible suppression of the beam instability by 
resonance broadening. The electron distribution f(v) is the sum of a thermal Maxwellian 
of temperature T e and a beam with velocity ty = 9 Vx e for Vx e the electron thermal 
velocity. The shaded grey, corresponds to the region of positive slope df/dv > 0 in the 
distribution function. It is centered around ty with a width A Vb- For a resonance width 
such that Av = v 2 — t’i < Au b , the growth rate 'y(Av) oc Af/Av > 0 and, hence, the 
instability proceeds, but at a slower rate as the effective slope is smaller. For an even 
larger resonance width, covering Av = V 4 — v :i > A Vb, we have j(Av) oc Af/Av < 0 and, 
hence, the instability is suppressed, i.e. the beam is stable to the generation of Langmuir 
waves. 
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